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A simple example of especially constructed potential barrier enables to show 
analytically (not numerically) the existence of tunneling effect for a Sommerfeld 
particle. 
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Let us remind that Sommerfeld model of charged rigid sphere [1] (Sommer- 
feld particle) is the simplest theoretical model to take into consideration the 
"back-reaction" of self-electromagnetic field on equation of motion of a radiat- 
ing extended charged body (for point-like charge we have the known Lorentz- 
Dirac equation with all its problems: renormalization of mass, preacceleration, 
runaway solutions, etc.). 

For the case of simplicity here we consider the nonrelativistic, linear in ve- 
locity, version of Sommerfeld model [2] . 

So let the total charge of a uniformly charged sphere be Q, mechanical mass 
- TO, radius - a. Then its equation of motion reads [1,2]: 

mij = Fe^t + V Wit - 2a/c) - v{t)] (1) 

here rj = -^^^ v = dR/dt, R - coordinate of the center of the shell, F^xt - 
some external force. 

This equation has no runaway solutions and solutions with preacceletation. 
But, as was shown in [3], it has the solutions, which one can interpret as classical 
tunneling. The physics of such effect is simple - Sommerfeld particle due to 
retardation begins to "feel" the action of potential barrier too late, when this 
barrier is overcome. 

In [3] the existence of solutions with tunneling effect was demonstrated with 
the help of numerical calculations. 

Here we present one simple problem, which enables us to see the appearance 
of classical tunneling without numerical calculations but using only analytical 
considerations. 

Let the Sommerfeld particle move along z -axis ( = (0, 0, i?) ) in a static 
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electric field E = (0, 0, Ez) produced by electric potential in a form of a step: 

(j) = B[e{z) - e{z - S)], Ez = = -B6{z) + BS{z - S) (2) 

az 

Then the external force Fext, acting on rigid sphere with density of charge 



Q 
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reads 



Fext = J dfpEz = ^ J dfS{\f- R\ - a)[-5{z) + 5{z - S)] 

With new variables f=£^+R, df = d^ = ^^d^ sin 9d9d(p, substituting cos 6 = ^ 
(it leads to z = R + ^n), integrating over ^ and ^, we get 
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j dn[-S{R+aii)+5{R+an-S)] 



QB 
2a 



R+a 



R+a-S 



dw5{w) 
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If S > a then the expression in square brackets yields the result 

OforR<-a, -1 for -a < R < +a, for +a < R < S-a, +1 for S- 

a<R<S + a, for R> S + a 

After this, eq.(l) in dimensionless variables y = -R/L, T = ct/L, d = 2a/L, 

for the relations S = L = 2a (taken for simplicity), is reduced to 



dT^ 

herefc=^, A=^, 
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dyij - 1) dyij) 



dr 



dr 
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and 



0, y<-l/2, 

-1, -l/2<y<l/2, 

+1, 1/2 < y < 3/2, 

I 0, 3/2 <y, 

Let us note that in the limit of zero particle's size a ^ the force F^xt tends 
to two delta- functions: 



lim Fext = QB[-5{R) + d{R - S)] 

a—^0 

This expression gives the following Newtonian eq. of motion for point-like par- 
ticle 

mv = QB[-S{R) + 6{R - S)] = -QB^[-eiR) + 9{R - S)] (4) 

uR 
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Equation ( 4) has the first integral: 



2 2 

my + QB[e{R) - e{R - S)] = const = 



Thus for < ii < 5 we have 



{2QB/m) 



This solution describes the overcoming of a point-like particle the potential 
barrier only if the initial velocity Vo is greater then the critical value Vcr, equal 
to 2QB/m or in dimensionless form 

^ = yer = \/2A (5) 

If fc = then the equation (3) for -1/2 < y < +1/2, < r < ii has the 
solution 

y = -1/2 + voT - Xt^/2, y = vo-XT (6) 

It follows from (6) that the point y = +1/2 is achieved (if the velocity y is 
always positive) at the moment of time ti : 



vo - ^{voy - 2X 
t, = ^ (7) 



Consequently the potential barrier is overcome if vq > v2A, and the minimal 
initial velocity for this is vo = V2X with y{ti) = and 

tl = vo/X = 2/vo (8) 

(the positive part of the force in the equation (3) for 1/2 < y < 3/2 does not 
hinder the motion of the particle, so for analysis of conditions for overcoming 
the barrier one can consider only the interval: —1/2 < y < +1/2 ), 

so for tl ~ 2, one must take vq = 1 and A = 1/2, and for ti = 3 — vq = 2/3 
and A = 2/9 and so on. 

If k in (3) is not equal to zero and if the point y = 1/2 is achieved at the 
moment of time t = ti in such a way that y{ti) = (this means that in the 
preceding moments of time the velocity is positive) then the further motion of 
particle, 1/2 < y < 3/2, is governed by the equation 



dr'' dr dr 

where / > 0. 

This eq. with the above condition y{ti) = has the solution 

r 

y = ef-'^-) / dt'f{t')e'-''''^ > 
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consequently the barrier is overcome and the particle continues its motion for 

T>tl. 

Taking all this into consideration, lets construct the " tunneling" solutions 
of the equation (3). 

Let the Sommerfeld particle reach the point y = —1/2 at the moment of 

time T = 0, J/ = 1/2 — at t = ti. 

Let the particle's velocity for r < be constant: 

dy 

— = vo = const 
dr 

Let the condition be true: y{ti) = (i.e. the barrier is overcome) and let 

tl = 2 (9) 

Then one must integrate the eq. (3) into two stages: in the first, for the 
interval 0<r<l, —l/2<y<yi and, in the second, for the interval 
1 < T < 2, yi < y < 2. 

Solution on the first interval is 

y = .o-^ + ^e(--) (10) 

and on the second — 

rC rC 

ao = ^(e('=") + l)-Ae('=^' (11) 
k 

Conditions y{ti = 2) = 1/2 y{ti = 2) = are reduced to the set 

Vo = /lA 
A = 1//2 

here 

f^ = -^+ p 2e(-^(l + i + ^) (13) 

The analysis of the functions fi{k), f2{k) (12, 13) shows that the condition 
Vo < V2X for the solution (10, 11) is vaHd due to the inequality /i / \/2/2 < 1 - i.e. 
there is the tunneling - the initial velocity has the forbidden, from newtonian 
point of view, values, but ought to the inequality /1//2 > 1 all of them are 
greater then 1: > 1 ( for /c = 0, in accordance with (8), Vo = I and A = 1/2 
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Let now instead of (9) be 

ti=3 (14) 

Then one must integrate the eq. (3) into tree stages: in the first, for the 

interval < r < 1, —1/2 < y < yi with the sohition in the form (10), in the 
second, for the interval 1 < r < 2, Vi < y < y2 with the solution in the form 
(11), and in the third, for the interval 2<t<3, 2/2<2/< 1/2 with the 
solution in the form 

y = v^-3^ + e(-'=*)(oi + ke^''\ao - Xe^'^hr + kXe^'^'^^T^ /2), 
k 

ai = e(2'=)(^ - 2A + 2kX - Xe^-''^ + ye^-''^) (15) 
k k 

Conditions y{ti = 3) = 1/2 and y{ti = 3) = are reduced to the set 

vo = hX 
X = 1/h 



here 



A-3k) 

-— — (6e(3*=) - {k^ + 2k + 2)e(2'=) - (4fc + 2)6^=) - 2) (16) 



p(-3fe) 

/4 = __((6fc+12)e(3'=)-(3A;3+7A;2+10fc+6)e(2'=)-(12fc2+i0A;+4)e('='-6fc-2) 

(17) 

Analysis of functions (16, 17) shows that the condition vq < V2X for the solution 
(15) is valid — fa/V^fi < 1 - i-e- there is the tunneling, and, contrary to the 
above case, the range of values of t^o is: 2/3 < t;o < 1 (for A; = in accordance 
with (8), = 2/3 and A ^ 2/9 ). 

The construction of tunneling solutions one can continue for ti = 4, 5, 6, ... 
and so on. It is obviously that the greater is the value of ti, the lower is the 
value of the minimal velocity vq necessary to overcome the barrier, but the more 
is complicated the form of the solution. 
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